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Abstract 

A completely quantum description of Raman process is used to investigate the nonclassical properties of the modes 
in the stimulated, spontaneous and partially spontaneous Raman process. Both coherent scattering (where all the 
initial modes are coherent) and chaotic scattering (where initial phonon mode is chaotic and all the other modes are 
coherent) are studied. Nonclassical character of Raman process is observed by means of intermodal entanglement, 
single mode and intermodal squeezing of vacuum fluctuations, sub-shot noise and wave variances. Joint photon-phonon 
number and integrated-intensity distributions are then used to illustrate the observed nonclassicalities. Conditional and 
difference number distributions are also provided to illustrate the nonclassical character. The mutual relation between 
the obtained nonclassicalities and their variations dependent on phases, rescalled time and ratio of coupling constants 
are also reported. 

Keywords: Raman scattering, nonclassicality, quasidistribution. 

1 Introduction 

Quantum statistics of Raman scattering were discussed from various points of view in a number of papers (see [TJ [2] 
and references therein, Section 10.4 of [3] and [I] for reviews). In this paper we describe the Raman scattering process 
with a completely quantum mechanical Hamiltonian. The model is capable to include the stimulated, spontaneous and 
partially spontaneous Raman process. We use the second order short-time approximation for solution of the Heisenberg 
equation of motion corresponding to this Hamiltonian. The solution is then used to relate the nonclassical properties 
of photons and phonons in these processes (i.e. stimulated, spontaneous and partially spontaneous Raman process). To 
be precise, nonclassical characteristics of photon and phonon modes generated in these processes are exhibited through 
a number of properties, including squeezing of vacuum fluctuations, quantum entanglement of modes, sub-shot noise 
and wave variances. Further, joint photon-phonon number, integrated-intensity distributions, conditional and difference 
number distributions are also found useful to illustrate the observed nonclassicalities. 

Remaining part of the paper is arranged as follows. In Section [2] we briefly describe the fully quantum Hamiltonian 
of Raman process and several criteria of nonclassicalities that are used in the present study. In Section [2] we have briefly 
described the model Hamiltonian that provides a completely quantum description of the Raman scattering process and 
also describe a normal-ordered characteristic function in the Gaussian form for the same system. Criteria for testing of 
nonclassicalities are also introduced. In Section [3] we investigate nonclassical character of Raman process for coherent 
scattering by means of intermodal entanglement, single mode and intermodal squeezing of vacuum fluctuations, sub-shot 
noise and wave variances. In Section|4]the same nonclassical characteristics are investigated for the chaotic scattering. The 
observed nonclassicalities are further illustrated through joint photon-phonon number and wave distribution in Section 
[5] In Section [6] we study difference and conditional number distributions associated with the Raman process. Finally, 
Section [7] is dedicated to conclusion. 

1 email: anirban.pathak@gmail.com 
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2 The model Hamiltonian and the criteria of nonclassicality 

A fully quantum description of the Raman scattering process can be provided by the following Hamiltonian [TJ [5] : 

H = hajjOjCij — (hgaLcJgcJy + hx* aL a vC L A + n - c -J > (1) 



j=L,S,A,V 



where h.c. stands for Hermitian conjugate and the subscripts L, S, A and V correspond to pump (laser), Stokes, anti-Stokes 



and vibration (phonon) modes respectively, Uj, dj and at are frequency, annihilation operator and creation operator in 
the j th mode, g and x are the Stokes and anti-Stokes coupling constants. Using the Hamiltonian ([I]) we can construct 
a set of Heisenberg equations and solve them in short-time approximation. A second order short-time approximated 
solution was already reported |3J. It is interesting to note that using the short-time approximated solution, we can obtain 
a normal-ordered characteristic function in the Gaussian form. Such a characteristic function can completely characterize 
the Raman process, and it can be analytically expressed as [5] 



+ (!q(t)$ + c.c.) 

+ E J<fc (D jk (t)P*p* k + b jk $)Pik + c.c.) }) 



2 

J I 

1 C* (t)/3] + c.c.) + frQ (t) - (i)] + (2) 



where c.c. stands for complex conjugate terms, the set (L, S, A, V) is assumed to be ordered and 



B L (t) = \x\ 2 t 2 \Z->'-' 



A 



B S (t) = i^m^i 2 , 

B A {t) = 0, 

B v (t) = \9\ 2 t 2 \t L \ 2 + \x\ 2 t 2 \ti 

C L (t) = -g*xt 2 tsUeM-1iu L t), 

C s (t) = 0, 

cu(t) = o, 

Cv(t) = -gxt 2 Z* s UeM-2iuvt), (3) 

D LS {t) = -±\g\ 2 t 2 £ L Z s exp[-i{uj L + uJs)t}, 

D LA (t) = -l\x\ 2 ^LUexp[-Hu L +u A )t] f 

D L v{t) = [ixt^A - \ (M 2 + |x| 2 ) i 2 khiv\ exp \—i {u L + u v ) t] , 

D SA {t) = -y X *t 2 &exi?[-i(u>s+UA)t], 

D sv (t) = {igfei - \\g\ 2 t 2 ^ v - gxfU&) exp [-i (u s + w) t] , 

D AV (t) = -k\x\ 2 t 2 U^vexp [-i {lo a + ui v )t] , 

D L s(t) = ~\gx*t £,U A exp [i (uj l - u> s ) t] , 

all other Dj k = 0; j — L,S,A,V are initial coherent complex amplitudes. As the above characteristic function is 
Gaussian consequently (3) can be used to obtain normal fluctuation quantities (variances) (^(AWj) 2 ^ and (AWjAW k ) at, 
which are defined as [3]: 

((AW-) 2 ) w = (af(t)a 2 (t))-( aj (t)a 3 (t)) 2 

= {B 2 + \C J \ 2 + 2B^(t)\ 2 +[C^f(t)+c.c.}) 

and 

(AWjAWh) N = (ot(t)ot(t)a i (t)ofc(t))-(at(t)o i (t))(oi(t)o fc (t)) (5) 
= (|^ jfc | 2 + \D ]k \ 2 + [D jk t*(t)C k (t) - D jk $j(t)e k (t) + c.c] ) . 

Brackets on the right-hand side in ([2] [4] and [5f mean an average over the initial amplitudes. Equations ([3])- ([s]) provide us 
with sufficient mathematical framework required for the study of the nonclassical character of stimulated and spontaneous 
Raman process. This is so because the criteria for various nonclassical phenomena can be conveniently expressed in terms 
of the quantities described in ([3])-([5]). For example, we may note the criteria for principle squeezing of vacuum fluctuations 
in single mode (j) and compound mode (ij), which are [3J: 

X j = l + 2(B j -\C j \)<l (6) 

and 

\ %j = 1 + B i + B j - 2ReD ij -\Ci + C 3 + 2£><j I < 1 (7) 

respectively. From the above two criteria it is clear that (|3| provides us sufficient input for analytic study of the principle 
squeezing of vacuum fluctuations, both in single modes and in compound modes. Similarly, condition for entanglement is 
in general [3] 
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{K lJ ) ± = (Bi ± \Ci\)(Bj ± \Cj\) - (| Ail T IA-I) 2 < 0, (8) 
and condition for sub-shot noise is 

Ch = B 2 + B) + IQI 2 + \C 3 \ 2 - 2|A,f - 2|Ai| 2 < 0; (9) 

further the condition for nonclassical sum- or difference-variance is 

((AW^-) 2 ) ± = ((A^) 2 )^ + ((A^-) 2 ) w ± (AW 3 AW k ) N < 0. (10) 

Present work aims to rigorously investigate the presence of different nonclassicalities in the Raman process in the second 
order short-time approximation. To begin with we will discuss intcrmodal entanglement in the next section. Before 
we present our analytic results it is important to note that for the convenience of understanding the process we have 
introduced following two scaled quantities: gt = t and W = p. The time evolution of various nonclassical characteristics 
can now be expressed with respect to dimensionless quantity gt — r and the ratio between the Stokes and anti-Stokes 
coupling constants p. Further we have used Ij = \£,j\ 2 for the incident stimulating intensities and the phases of the complex 
amplitude £j = \ exp{(f>j) are denoted as <f>j that are combined as 

4>L - 4>v - 4>s = #2 

and 

§A - <j>L - <t>V = #1 ■ 

where #2 and 9\ can be visualized as the mismatch phases in Stokes {us = wl — wy) and in anti-Stokes {ua = wl + uy) 
transitions, respectively. In the following the coupling constants g and \ are assumed to be real. 



3 Phonon mode is coherent 

In the above discussion, all the modes including the phonon mode are coherent. In such a situation we may investigate 
the existence of different kind of nonclassicalities by using (|3])-(|To|) . The same is done in the following subsections, where 
intermode entanglement, single mode and intermode squeezing, sub-shot noise and variances are studied in detail. 



3.1 Intermodal entanglement 

Substituting ^ in ^ we obtain 

(K L v)+ = (K LV )_ = -p 2 r 2 I A = E LV (11) 

and 

(K sv ) + = (K sv )_ = -t 2 I l = E sv . (12) 

(Kij) ± — for all the other cases. Since we are using a second order short-time approximated solution we cannot conclude 
anything about the separability of those four modes for which (Kij) ± = 0. But we can clearly see that in stimulated 
Raman process (where I a 0, II =/= , Is =/= 0, ly ^= ) the vibration-phonon mode is entangled with the pump-mode 
and the Stokes mode and it does not depend on Is and ly. Consequently if we think of a partially spontaneous Raman 
process with I a ^ 0, II 0, Is = 0, ly = 0, then also we will observe both type of photon-phonon entanglement that we 
have observed in stimulated Raman process. Interestingly in the completely spontaneous process (where I a = 0, II ^ 0, 
Is = 0, ly = 0) we can also observe entanglement between Stokes mode and phonon mode, but in such situation we 
cannot conclude about the separability of the pump mode and the phonon mode. 



3.2 Single mode and intermodal squeezing 

Substituting Q in ( ^ and Q we obtain in the interaction picture reflecting the compensation of exponential function 
in (pi in a homodyne detection 
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Figure 1: Variation of — 1 with respect to p and gt when = 10, \£a\ = 1, = 9, = 0.01 and <pL — 4>A- 



From the above equations one can easily obtain the following conditions: 

1. Since > \^a\ in general the pump mode is always squeezed if |<7| > otherwise it is squeezed if p\£,a\ < |£s| 5 
which is expected to be satisfied in most cases. 

2. \la < 1 in stimulated Raman process if \^a\ < + |£s|), which is the natural case. So intermodal squeezing 
between pump and anti-Stokes mode can be observed in stimulated Raman process. In spontaneous Raman process 
Xla — 1 so squeezing is not observed, but in partial spontaneous Raman process with = 0, |£s| = 0,0 < < 

squeezing can be observed. 

3- Xsa < 1 iff p > 1, i.e. if anti-Stokes coupling is stronger than the Stokes coupling. If p > 1 then the intermodal 
squeezing in Stokes and anti-Stokes modes is observed for both stimulated and spontaneous Raman processes. 

4. For a completely spontaneous process Xav ~ 1 + t 2 Il is always greater than 1. However, also in the stimulated 
process the term t 2 Il will be dominant. The same is the case for Xls- 

5. For a very short time the linear term in X^y would dominate and consequently, during that time Xlv ~ 1 — 2prlA 
will be less than unity and consequently squeezing will be observed in stimulated and partially spontaneous process. 

6. For a very short time the linear term will dominate in Xsv an d consequently, during that time Xsv ~ 1 — 2t|^| < 1 
would indicate intermodal squeezing in both stimulated and spontaneous process. 

Variation of — 1 with respect to p and r are shown in Fig. [T] which clearly depicts the above observations. 



3.3 Sub-shot noise 

When we substitute (|3| in |)| we find that 

C LV = ~2p 2 T 2 I A = 2(K LV ) ± (14) 



and 

C S v = -2t 2 I a = 2{K SV )± (15) 

and Cij = in the remaining four cases. Since in this particular system C^v and C$v are directly proportional to (Kly)± 
and (Kgy)±, wherever we have seen intermodal entanglement we, can also observe sub-shot noise in those cases. 



3.4 Variances 

Using pl-dsl and (10 1 we can obtain the analytic expressions for intermodal variances in the following forms 



(Aw Av f) = 2T 2 i v (i L ~ P \^s\\U\cos{e 1 + e 2 )), 

(AWsa) 2 ) = 2T 2 I L (I s -p\i s \\U\cos(8 2 -6 1 )), 
( + 

(AW LA ) 2 ) = -2 P T 2 I L \is\\U\cos{9 2 -6 1 ), 
( + 

(AW LS f) + = 2pr 2 I L \U\(p\U\ + 008(02-0!)), (16 ) 

(AW LV f) = 2 p 2 t 2 (M l I a + M a I v + Ia-IlIv) 

+ 2 P T 2 \Z S \\U\ (II cos(0 2 - 0i) + I v 008(0! + 2 )) + 4|>t|£a||&||£v| cos(0i), 

slv) 

2pT 2 \ds\\U\ (21 l cos(0 2 - 0i ) - 3/ y 008(0! + 2 )) + 4t|&| \Zv\ cos(6 2 ), 



(AWsv) 2 ) = 2t 2 (3I l I s + 3I l I v +P 2 IaIv + Il- Isl- 



and 



{AW Av f)_ = 2T 2 I v (l L + 2p 2 I A -p\^ s \\U\cos(9 1 +e 2 )), 

(AW SA ) 2 )_ = 2T 2 I L {I s +p\(is\\U\cos{e 2 -e 1 )), 

(AW LA f)_ = 2 P T 2 I L (2pI A ~ |esll^|cos(0 2 -0i)), 

{AW LS f)_ = 2T 2 I L ( P 2 I A + 2I s -3p\Zs\\U\cos(e 2 -ei)), (17) 

(AW LV ) 2 )_ = -2p 2 r 2 (I L I A + I A I V +I A - I L I V ) + At 2 I l I v 

- 6pr 2 \^s\\U\ (h cos(0 2 - 00 + I V cos(0! + 2 )) - A P t\U\\U\Cv\ 008(9^, 
(AW SV ) 2 )_ = -2t 2 (I l I s + I l I v -p 2 I a I v +I l -I s Iv) 

- 2pT 2 \^s\\U\ (21 l cos(0 2 -Oi)-I v 008(0! + 2 )) - 4r|a||es||Cv| cos(0 2 ). 

Negativity of intermodal variances ((AW)- 3 ) ± implies nonclassicality. Analytic expressions for intermodal variances 
((AW) 2 j) + and ((AW)^-)_ for all the possible combinations of modes in the stimulated Raman process are provided in |l6| 
and (|17|) , respectively. It is difficult to conclude directly about the presence of nonclassicality from these general analytic 



expressions of ((AW)|.-) + . Thus to visualize the existence of nonclassicality we have plotted the analytic expressions 
provided in (16 1 and ( jl7| . The plots are given in Fig. |2j and Fig. |3j and it is easy to see that both ((AW) 2 j} + and 
((AW) 2 j)_ depicts nonclassical behavior for a) pump and phonon mode and b) pump and anti-Stokes mode, c) Stokes 
and phonon mode. However for pump and Stokes mode only ^(AT / F)| /S ) + shows the existence of nonclassicality. 

For the chosen values of and p we have seen that ((Alf)j A V = 2t 2 Il\£,s\ (|£s| ± p\£a\ cos(0 2 — 0i)) > but a 
negative value is possible if p\£ A \ > \£,s\- Thus a very strong anti-Stokes coupling (compared to Stokes coupling) may yield 
nonclassical variance for Stokes and anti-Stokes mode. This is consistent with the appearance of intermodal squeezing 
where A as — 1 = Ilt 2 (1 — p) is negative only when p > 1, that is when anti-Stokes coupling is stronger than Stokes 
coupling. 

Now from |l7| ) we can easily observe that for a completely spontaneous Raman process ^(APF)| V )_ = — 2t 2 I[ j is 
always negative which indicates intermodal nonclassical behavior between phonon mode and Stokes mode. We have 
already shown that these two modes show intermodal entanglement, sub-shot noise behavior and squeezing of vacuum 
fluctuations in the spontaneous Raman process. Thus as far as the nonclassicalities in spontaneous Raman process are 
concerned these two modes play the most important role. 



From p6[ ) and ( 17 1 we can see that for very small values of rescaled time r the term linear in r is expected to dominate 
in ((AW)]~^) ± and in ((AT / F)| y ) ± ; ((AW) 2 LV ) ± varies with 0i, which is exhibited in Fig. |ij Further the linear term in 
((AVF);|y) + is very weak and the nonclassical behavior can be seen only for a very small values of r. This is why in Fig. 
[2] we have plotted {(AW) 2 SV ) + for a very short time only. 

4 Phonon mode is chaotic 

In this case we perform the average over the initial phonon amplitude in ^ with a Gaussian distribution in the Gaussian 
approximation. Assuming that the phonon mode is chaotic with average phonon number (ny) , then the coefficients in 



Figure 2: Variation of intermodal variance ((Alf)?-) with respect to gt and 02 here we have chosen — 10, = 1, 
N = 9, 1^1=0.01, and p= 0.9. 




Figure 3: Variation of variance ((Alfjj^ with respect to gt and 02 here we have chosen = 10, = 1, \£s\ = 9, 
= 0.01, 6>i = |, and p = 0.9. 




Figure 4: Variation of variance ^(AT / F)| / y) ± with respect to 0! when Q 2 = \ a U other parameters are same as in previous 
figures. 



the interaction picture described in |3| get modified as 



B L 


= \x\ 2 t 2 \u\ 2 i(n v ) + r 

= \9\ 2 t 2 \U 2 ((nv) + I) 

= lxl 2 * 2 l^l 2 K-), 


+ \g\ 2 t 2 \ts\ 2 (nv 


B s 




B A 




By 


~ ( n v), 




C L 


= -9*Xt 2 UZs(2(nv)^ 


-1), 


C s 


= o, 




C A 


- o, 




C v 






D LS 


= -\\g\ 2 t 2 ^Us(2(n v ) 
= -hx\ 2 t 2 aU(2(n v ) 
= -bx*t 2 &(2(n v ) + 


+ 1), 


D L a 


+ 1), 


Dsa 


1), 


D S v 


= igt£,L ((nv) + 1) j 




D L v 


= ((nv) + 1) ! 




Dav 


= o, 




Dls 


= -9X*t 2 e A (L ((nv) + 


1), 


D L a 


= -9X*t 2 (isZL(nv), 




Dsa 


= o, 




Dsv 


= o, 




D L v 


= igt£ s (n v ), 




D av 


= ixtZl(n v )- 





(18) 



Now using equations Q , (|5| and ( 18 1 we can obtain for single-mode variances 

((AW L f) N = 2p 2 r 2 I A I L ((n v ) + 1) + 2T 2 I s I L (n v ) - 2 P T 2 \U\\ts\h (2(n v ) + 1) cos(# 2 - 6 1 ), 

((AW s ) 2 )n = 2T 2 I L I S ((nv) + l), riql 

((A^) 2 )^ = 2p 2 r 2 I L lA(nv), { ' 

((AW v f) N « (n v ) 2 , 

and for correlation fluctuations 

(AW L AW A ) N = -T 2 P 2 I L lA(2(nv) + l) + 2pT 2 lL\is\\U\(nv)cos(e 2 -e 1 ), 

(AW L AW S ) N = -T 2 lLls(2(nv) + l) + 2 P T 2 I L \is\\U\((nv) + 1) cos(0 2 - (20) 

(AW S AW A ) N = -T*pI L \U\\ts\ cos(8 2 - 9 1 ) (2(n v ) + 1) . 

Analytic expressions of the other cross-correlations are not of interest as all variances that involve phonon mode will 
always be positive because of the dominance of ((AW-V) 2 ) « (nv) 2 term. Now substituting equations (18 1, (19 1 and 



(20 1 in the criteria of nonclassicalities introduced in (|6|-([10j) we can investigate the nonclassical character of stimulated 
and spontaneous Raman process when the phonon mode is chaotic and then compare the results with the similar results 
obtained in the coherent case. This is done in the following subsections. 

4.1 Intermodal entanglement 

Substituting (18 1 in ^ we obtain 

(K SV ) ± = -r 2 I L ((n v ) + l), , . 

(K LV ) ± = -p 2 T 2 I A ((nv) + l)T2 P T 2 \Ums\(nv)(5(n v )+2). ^ 

All other (K ij ) ± = 0. 
We can conclude: 



1. From (21 1 it is clear that the phonon mode is always entangled with Stokes mode. The same characteristic was also 
observed in coherent case but if we consider (Ksv) + 35 a- measure of amount of entanglement, then the amount of 
entanglement in chaotic case is increased by a factor of (1 + (nv)) and it is more announced. 

2. Similarly the phonon mode can be entangled with the pump mode. It is straightforward to see that (Klv) + = 
—p 2 r 2 I A ((ny) + 1) — 2pr 2 |£ J 4||£s|(ny) ((3(ny) + 2)) exhibits intermodal entanglement. 

3. But interestingly (Klv)_ does not show signature of intermodal entanglement. 

4. Stokes mode and phonon mode are entangled for completely spontaneous Raman process also but the present 
calculation is non-conclusive about entanglement of pump and phonon mode. 
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4.2 Single mode and intermodal squeezing 

By substituting Q in M and we obtain 

Al = l + 2T 2 [p 2 I A ((n v )+l)+I s (n v )-p\UMs\(2(nv) + l)], 



X LA « l+p 2 T 2 |a|(|ai(«y)-|e4|(2(m/) + l)), 
A iS « l + r 2 |a|(|eL|((ny) + l)-|esl(2W> + l)), 
Asa = l + r 2 / L (l-p)((n v ) + l-p<ny)). 



We see that: 



(22) 



1. Squeezing in the pump laser mode can be observed by approximating > ( 2 {n v )+i) ■ 

2. Intermodal squeezing is not possible when one of the mode is phonon mode as in that case A^y rj 1 + (nv) > 1- 

3. Intermodal squeezing will not be usually observed between pump mode and anti-Stokes mode as Xla < 1 implies 

< 2|£a| + ~ But technically it is allowed and intermodal squeezing between pump mode and 

anti-Stokes mode can in principle be seen for stimulated Raman process as well as for partially spontaneous Raman 
process (I A ? 0, I L ? 0, I s = 0, I v = 0). 



4. Intermodal squeezing between pump mode and Stokes mode is possible if ((ny) + 1) < \£s\ (2(ny) + 1) , i.e. if 
iff! < ?(»v ) +1? • For ( nv ^ ^ 1 this condition implies that |£ L | < 2|£g| and for (ny) = it implies |£ L | < |£g|- 

5. To have \$a < 1, we need p = l# > 1 and ((ny) + 1) > p(ny), i.e. W < 1 + t^t- For coherent scattering (ny) — > 
and we have the previous condition M > 1. 

6. Similarly when p < 1 then Asa < 1 implies (p(ny) — ((ny) + l))>0=>p>l + => p > 1. Thus the condition 
of negativity is not satisfied and we are non-conclusive about the entanglement between Stokes mode and anti-Stokes 
mode if p < 1 . 

4.3 Sub-shot noise 

By substituting ([ll) in ^ we obtain 

Cav = (ny) 2 (l-2p 2 r 2 / L ), 

C LV = (n v ) 2 -2p 2 r 2 I A {(n v ) + l) 2 -2r 2 I s (n v ) 2 , (23) 

Csv = (ny) 2 -2r 2 / L ((ny) + l) 2 , 

and all other dj — 0. For stimulated Raman process, sub-shot noise is observed in the above three cases. In coherent 
case subshot noise behavior was not observed for anti-Stokes and phonon mode. Further, negativity of Cav & n d Csv will 
be observed for spontaneous Raman process too. But in the spontaneous Raman process sub-shot noise behavior will not 
be observed for pump and phonon modes. However, we can observe it for partially spontaneous process ((ny) ^ 0, 4 ^ 
0, I A = and I s = 0.) 

4.4 Variances 



By substituting ( |19| and (|20| in (|10| we obtain 

) sa 



((AW) SA ) ± = 2p 2 T 2 I L I A (n v )+2r 2 I L I s ((n v ) + l) 

T 2T 2 pI L \U\\^s\ cos(9 2 - 9,) (2(ny) + 1) , 

((AW) 2 LS ) ± - 2 P 2 t 2 I a Il ((ny) + 1) + 2r 2 / L /s (2(ny) + 1) (1 T 1) m] 

- 2pT 2 I L \Z s \\U\cos(6 2 -e 1 )((2(n v ) + l)T2((nv) + l)), 1 ' 
((AW) 2 LA ) ± = 2p 2 T 2 / A / L ((2(ny) + l)T((ny) + l)) + 2r 2 /s/L(ny) 

- 2 P t 2 \U\MIl cos(6 2 - di) ((2(ny) + 1) =F 2(n v )) ■ 



From (24 1 we observe following: 

1. ((AW) LS ) + = 2pt 2 \£a\Il (p\£a\ ((ny) + 1) + |£s| cos(6> 2 - #i)). Thus negative variance can be seen forp|£^| ((ny) + 1) < 
|£s|- Since < |£s| this nonclassical feature between pump mode and Stokes mode will be observed for small 
values of mean phonon number (ny). 

2. In the analytic expression of ((AW) AS )± if we assume Ia <C Is then we obtain 

((AW) 2 AS ) ± = 2r 2 I L \^ s \ (l&l ((nv) + 1)TP\U\ (2{n v ) + 1) cos(9 2 - 9,)) 

which would show nonclassicality if Pj^j > 2 ™nv)+i • This implies p\£a\ > §|£s| which is inconsistent with the 
assumption I a <C Is- Thus if I a <C Is then we do not observe nonclassical variance in anti-Stokes and Stokes 
modes. 



Figure 5: a) Joint photon-phonon number distribution for Stokes and vibration modes (left) and b) Quasi-distribution 
of integrated-intensities for the same modes (right). Here Bg = 0.1, Bl = 0.01, By — 0.11 and s = 0.7. 



3. If we assume I a ^ Is an d consider the complete analytic expression of {{AWas) 2 )±, then the condition > 
2{nv)+i wm serve as necessary but not sufficient condition of nonclassicality. Now if we assume that for some choice 
of p, I a. Is? ( n v) we observe nonclassical intermodal variance for Stokes and anti-Stokes modes, then we can show 
that for that situation {{AWls) 2 ) + will not show nonclassicality. The proof is simple. First we assume that both 
((AW LS ) 2 )+ and ((AW AS ) 2 )± are negative. Therefore, j^y+i > p]ff( > afcvH-i ' which implies ( 2 (™v) + 1) > 
{{ny) + I) 2 or {ny) 2 < 0. Thus by reductio ad absurdum we have shown that intermodal nonclassical variance 
cannot be seen simultaneously in a) Stokes and anti-Stokes mode and b) Stokes and pump mode. 

4. For the compound mode {LA) one could observe sub-shot noise provided that |xI|£a| > l/l • 

5 Joint photon-phonon number and wave distribution 

We can illustrate the above results for nonclassical behavior of modes in Raman scattering by joint photon-phonon number 
and integrated-intensity distributions along the lines given in [6 (and references therein) in Gaussian approximation. 
For simplicity we consider scattering by phonon vacuum (in optical region and for room temperature {ny) w 0) for 
compound modes {SV ) and {LV), which exhibit quantum entanglement up to the second order in t. From |3| we see 
that Kgy = —Bg = —By = — \g\ 2 t 2 lL provided that we consider spontaneous scattering {Ia = Is = 0; in this case 
Cg — Cy = Dgy =0). In principle we can also consider partially stimulated scattering with Ig 0, when using shifted 
distributions in Wg along Ig [7] to adopt spontaneous process. Thus Kgy + Bg = Kgy + By = and from the formulae 
given in [5] we obtain the joint photon-phonon number distribution 



{Bg) ns 

p{ng,n v ) = ————— 5 ns , nv , (25) 



(l + B s ) 1+ns 

i.e. it is diagonal expressing a pairwise structure of photon-phonon process in this case. It is shown in Fig. [5]i. The 
corresponding s-order quasidistribution of integrated-intensities is [5] 

sinf w r Wv 

P s {Wg,W v ) = e-^T- y • ' , (26) 

nBg s Wg - Wy 

where Bg s = Bg + an d Kgy s = Kgy + (1 — s)Bg + s is ordering parameter. For the threshold value of the 

ordering parameter we have s t h = 1 + Bg + By — {Bg + By) 2 — 4Kgy sa 1 + 2Bg — 2\/Bg. Choosing \g\t = 0.1, II = 10, 
we have s t h — 0.57. So we calculate the quasi-distribution for s — 0.7; in this case Bg = 0.1, Bg s = 0.25, Kyg\ s = —0.048. 
This quasi-distribution is shown in Fig. J5|). It takes on negative values exhibiting nonclassical oscillations and behavior. 

Similarly, we can treat the compound mode {LV) considering again photon vacuum scattering with partial stimulation 
Ia 7^ and Ig = 0. Shifting distribution in Wl along 1^ and neglecting short-time terms as above we obtain from 
(j), K LV = z \x?t 2 I A = -B L and By = \g\ 2 t 2 I L + \ X \ 2 t 2 I Al i.e. K LV + B L = and K LV + By = \g\ 2 t 2 I L > 
{Cl = Cy = Dlv = 0). Thus for the joint photon-phonon number distribution [B] we obtain 

, s nyl {By + K Lv r-- n - 

p{n L ,n v ) = — - — 7—^ ,n v >n L . (27) 

n^-yny - n L )\ (1 + Bn v ) nv+L 

For ny < the distribution is zero. Its quantum behavior is illustrated in Fig. [6^, showing one-side behavior along the 
diagonal compared to the earlier cases [7]. For the threshold values of the ordering parameter we have s t h ~ 1 + \g\ 2 t 2 lL — 
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Figure 6: a) Joint photon-phonon number distribution for pump and vibration modes (left) and b) Glauber-Sudarshan 
quasi-distribution of integrated-intensities for the same modes with Wy — 0, 0.1, 0.2. • • • 2 (right). 




Figure 7: Conditional Fano factor Fy,c f°r phonon mode shows nonclassical behavior as Fy t c < 1- Here By = 0.11 and 
Sz, =0.01, K LV = -0.01. 



2\x\t\/lA- Assuming for simplicity \x\ = \g\, Ia = 1 and II = 10, we have sth = 0.9 and nonclassical behavior of wave 
quasi-distribution is illustrated by the Glauber-Sudarshan quasi-distribution of integrated-intensities for s = 1: 



sm 



Pn(W l ,W v ) 



1 



Wl »'v 

'IB,, IB, 



V 



-VF V 



(28) 



as shown in Fig. [6b. The existence of nonclassical character is clearly visible through the negative values of Pn(Wl, Wy)- 



6 Difference and conditional number distributions 

In this section we can further illustrate the observed nonclassicalities via difference and conditional number distributions. 
For example, nonclassical character associated with a mode can be illustrated using conditional Fano factor, which is 
defined as 

<(An,) 2 ) c 



F 



i,C 



for mode i. Corresponding condition for nonclassicality is F^c < 1. Analytic expressions for conditional Fano factor are 
obtained here for modes of interest (i.e. for Fl : c an d Fv,c) as follows: 



F, 



L,C 



1 



Bl 
By 



and 



F 



v.c 



(n L + 1) ( 



1+Bz 



- 1 



(n L + 1) 



l + By 

1+Bi 



(29) 



(30) 



It is now easy to observe from (j3| that Bl and By are always positive, consequently the conditional Fano factor Fl,c is 
always less than unity. Thus conditional Fano factor always depicts nonclassicality in pump mode. However, in phonon 
mode the presence of nonclassicality (i.e. Fyc < 1) i s n °t directly visible from the expression, but the same is shown in 
the Fig. [7] 



in 



Figure 8: Conditional number distributions pc (n^ny) (left) and pc {nv;nL} (right). Here By = 0.11 and Bl = 0.01, 
K LV = ~0M. 




Figure 9: Difference number distribution p_ (n) and p po i S (n). Subpoissonian character is shown by p-(n). Here By = 0.11 
and B L = 0.01, K LV = -0.01. 



The corresponding number distributions are obtained as 



Pc(n L ;n v ) = (l - (g^P , 

Pc{n v ,n L ) - nL <( nv _ nL y 1+By y i +Bv J \B V -B L ) 
These conditional number distributions are plotted in the Fig. [8] Difference number distribution can be obtained as 



(31) 



(By - B L ) U 

p - {n) = { i +Bv -B L r^ (32) 

((An)V - (B v - B L ) (1 + By - B L ) 



and Poissonian distribution for the same two modes is 



PPoUn) = ( Bv+ i B ^ n e-(*v + B L) . (33) 



A joint plot of P-(n) and p po i S (n) is provided in Fig. [9j which clearly shows subpoissonian character in p_(n). Thus 
a nonclassical difference number distribution is observed. For the sub-shot noise parameter R — ,^ A J\' 3 ? Is we have 
Ral- „ 2 % = 0.83 < 1. 



7 Conclusion 

We have observed different type of nonclassicalities in the stimulated, completely spontaneous and partially spontaneous 
Raman process. The observations that are discussed in detail in Section |3.1| are summarized in Table [T] for coherent 
scattering. We see that in general various nonclassical features of the process can or cannot be directly related, only 
for combined modes (LV) and (SV) all of them occurs simultaneously. We have not restricted ourselves to the study of 
coherent scattering alone. In Section [4] we have investigated various nonclassical characters of Raman process when the 
phonon mode is chaotic. Finally we have illustrated our results by joint photon-phonon number and wave distributions. 



1 1 



Mode 


Eij 




Cij 


/(AWij) 2 \ 

\ 1 — 


((AW-) 2 ) 




{Kij)- 


AV 


non- 
conclusivt 


>1 


non- 
conclusive 


+ve 


+ve 


non- 
conclusive 


non- 
conclusive 


AL 


non- 
conclusivt 


<1 if 
Il>Ia 
(expected) 


non- 
conclusive 


-ve region 
exists 


-ve region 
exists 


non- 
conclusive 


non- 
conclusive 


AS 


non- 
conclusive 


<1 if 

\x\ > M 

1 /V 1 1 J \ 


non- 
conclusive 


+ve 


+ve 


non- 
conclusive 


non- 
conclusive 


LS 


non- 
conclusivt 


>1 


non- 
conclusive 


+ve 


-ve region 
exists 


non- 
conclusive 


non- 
conclusive 


LV 


always 
-ve 


< for short 
time 


always 
-ve 


-ve region 
exists 


-ve region 
exists 


always 
-ve 


always 
-ve 


SV 


always 
-ve 


< for short 
time 


always 
-ve 


-ve 


-ve 


always 
-ve 


always 
-ve 



Table 1: Negativity of different characteristics of nonclassicality. Here we have used = BiBj — \Dij\ 2 1 A^- = 1 + + 
Bj - 2ReD ij - \d + Cj + 2£> y | < 1, = B? + B] + \Q\ 2 + \Cj\ 2 - 2|D y | 2 - 2|D y | 2 < 0, (K l0 ) ± = {B t ± ^{B, ± 

|C<|)-(|A,|T|A,|) 2 . 
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